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ABSTRACT 

The  high  pressure  supercritical  facility  at  AFRL  is  used  to  investigate  effects  of  chamber  pressure ' 
from  the  thermodynamic  subcritical  to  supercritical  values  at  a  supercritical  chamber  temperature.  At  subcntical 
pressures  the  jets^exhibit  wave-like  structures  which  amplify  downstream  and  eventually  break  up  ^t0 
shaped  small  entities.  Tire  formation  of  many  droplets  is  seen  at  higher  pressures,  resembling  a  second  wind-induced 
liquid-jet  breakup  Further  increase  of  chamber  pressure,  near  the  critical  condition,  fails  to  induce  the  transition  in  o 
foil  liquid  atomization  regime.  At  this  point,  the  jet  anatomy  changes  abruptly  to  mutate  turbulent  gas  jet  mjectiom  The 
to  SSZh  rate  is  plotted  against  the  chamber-to-injectant  density  ratio,  along  with  available  date  on  other 
liquid/gaseous  jets  and  mixing  layers,  producing  a  unique  and  informative  graph.  For  suPer“l  m  ?^d  C  OseW 
measured  growth  rate  agrees  well  with  a  theoretical  equation  proposed  by  Papamoschou  and  Roshko  [1]  and  c  os  y 
follows  thf  trend  of  Dimotakis  [2]  for  incompressible  but  variable-density  gaseous  turbulent  nuxmg  layers^  Fra 
isis  of  Set  interface  also  kU  a  similarity  to  gas  jet  behavior  ™th  comparable  “  “  £ 

first  time  quantitative  evidence  has  been  provided  to  support  qualitative  visualizations  suggesting  that  supercritical  jets 

appear  to  behave  like  conventional  gas  jets. 

INTRODUCTION 

A  general  trend  to  operate  under  an  increasingly  higher  combustor  pressures  is  observed  in  rockets  gas  terbmes,  and 
diesel  engines  primarily  due  to  enhanced  effects  on  thrust,  power,  or  efficiency.  Under  such  conditions,  the 
environment  into  which  fuel  is  injected  can  be  at  a  supercritical  state.  This  motivated  a 

AFRL  to  initially  document,  analyze,  and  finally  understand  anatomical  changes  of  a  liquid  cryugemc  jet  subjected  to  a 
transition  from  subcritical  to  supercritical  conditions,  as  occurs  m  the  thrust  chamber  environment  DU? 

chemical  rocket.  As  examples,  the  Space  Shuttle  main  engine  thrust  chamber  pre|ure  is  about  ^3  MPa^ 
combustion  chamber  pressure  for  Vulcain  (Ariane  5)  with  liquid  H2  (P0=1.27^a,  Tc~33.2^)tiiqui ld 
5  043kpa  Tc=154.58k)  can  reach  up  to  10  MPa  while  a  record  pressure  of  nearly  28.2  MPa  has  been  reported.  Very 
little  information  is  available  on  liquid  jets  injected  into  supercritical  environment. 

For  a  pure  substance,  the  distinction  between  liquid  and  gas  phases  disappears  at  and  above  the  critical  P^sure  (cailed 
as  “fluid”)  and  density,  thermal  conductivity,  mass  diffusivity,  and  surface  tension  show  large  variations  near  th 
critical  point.  At  elevated  ambient  pressure,  the  solubility  of  gases  into  the  fctiid  phase  increases  and  mulricomponent 
phase  equilibrium  informationf  or  “critical  mixing  temperature  or  pressure}  should  be :  used, ^see  Bruno 
Hence  when  a  pure  hquid  fuel  drop  is  introduced  into  a  gas,  a  thin  layer  that  is  a  mixture  of  that  gas  and  %  ^el  'S 
formed  on  the  Lace  that  spreads  spatial^n  time.  In  what  follows,  the  terms  submitted and 
reduced  pressure  (Pr)  and  temperature  CUfare  referenced  with  respect  to  the  critical  condition  of  the  pure  substance 
used  in  the  drops  or  jets,  and  not  the  environment. 

Recently,  some  results  on  liquid  jet  injection  into  supercritical  condition  have  been  presented,  for  example  Newman  rnd 
Brzustowski  [4]  and  Mayer  el  al.[5].  [6],  mi  Chen  mi  Sul  [7]  *  high  Reynolds  numbers  and  IRoodW mi  TaUey^ 
[81  at  low  Reynolds  numbers.  Newman  and  Brzustowski  [9]  proposed  the  possibility  of  gas 
supercritical  temperatures  and  pressures  the  jet  can  be  considered  as  a  variable-density  single-phase  turbukn 
submerged  gas  jit.  Also,  assuming  self-preserving  flow,  negligible  gravity,  zero  latent  heat  of  vaporization,  ideal  gas 


distribution  statemen 

Approved  for  Public  Releas 
Distribution  Unlimited 


behavior,  and  thermal  equilibrium  between  gas  and  drops,  they  developed  a  model  for  predicting  the  profile  of  the  outer 
extent  of  a  supercritical  steady  jet.  Comparison  with  experiment  was  poor  near  the  injector.  Mayer  et  al  [6]  used  a 
liquid  N2  (LN2)  jet  at  105  K  into  a  N2  environment  at  300  K  but  at  varying  ambient  sub-  to  supercritical  pressures  and 
observe  that  the  jet  behaves  similar  to  the  classical  atomization  of  liquid  fuel  with  ligaments  and  drops  below  the  critical 
pressure.  They  report  gas-jet  like  appearance  with  no  evidence  of  droplets  at  supercritical  conditions.  For  an  N2-into-N2 
cryogenic  jet,  Woodward  and  Talley  [8]  showed  that  the  addition  of  He  gas  to  the  supercritical  ambient  N2  fluid  forced 
the  jet  to  look  like  injection  into  the  subcritical  ambient  condition. 

EXPERIMENTAL  SETUP 

The  stainless  steel  chamber  in  Fig.  1  can  withstand  pressures  and  temperatures  of  up  to  137  atm  (2000  psi)  and  473  K, 
respectively,  and  has  two  facing  circular  sapphire  windows  and  two  UV-grade  side-mounted  slot-shaped  quartz 
windows  for  laser  in/out  of  the  chamber.  Transport  of  the  liquid  through  the  lines  and  the  injector  is  achieved  as  shown 
in  Figure  1.  The  liquid  N2  from  a  dewar  is  used  to  cool 
and/or  liquefy  the  injectant  passing  through  the 
cryogenic  cooler  prior  to  injection.  In  addition  to  N2,  a 
branch  is  connected  for  possibly  introducing  other  gases 
or  mixtures  of  gases  as  a  chamber  medium.  The  mass 
flow  rate  of  the  injectant  is  regulated  and  also  measured 
via  a  mass  flowmeter,  and  a  precision  micrometer  valve. 

For  more  details  refer  to  Woodward  and  Talley  [8]  and 
Chehroudi  et  al.  [9],  The  injection  in  this  study  is 
through  a  sharp-edged  5Q-mm  long  stainless  steel  tube 
with  1.59-mm  (l/l^”'8i’ameter  and  a  254-micron 
(0.0 1C£)  "inner  hole  (length-to-diameter  ratio  of  200). 

With  the  Reynolds  number  ranging  from  25,000  to 
70,000,  the  entrance  length  from  50  to  100  is  needed,  see 
Schlichting  [ 10] .  Tire  length  is  therefore  long  enough  to 
ensure  fully-developed  turbulent  pipe  flow  at  tire  exit. 

The  rig  is  fully  instrumented  with  thermocouples, 
pressure  gages  and  transducers,  and  mass  flowmeters  at 
tire  locations  indicated  in  Fig.  1.  Back-illumination  and  a 
model  K2  Infinity  long-distance  microscope  is  used 
with  a  TM-745E  high  resolution  (768(H)  x  493(V) 
pixels  in  8.8(H)x6.6(V)  mm  actual  sensing  area) 
interlaced  CCD  camera  by  PULNix  to  form  images  of 
tire  injected  jets.  For  more  details  refer  to  Chehroudi  et 
al.  [9J. 

ANATOMICAL  CHANGES  OF  THE  JET 

Figure  3  shows  images  of  the  N2  jet  injected  into  N2  at  a  fixed  supercritical  chamber  temperature  but  varying  sub-  to 
supercritical  pressure.  At  the  lowest  subcritical  chamber  pressure,  the  jet  is  liquid-like  with  surface  instabilities  that 
grow  downstream.  Very  fine  drops  are  seen  ejecting  from  the  jet  at  Pr  of  0.83  and  the  jet  grows  away  from  the  injector. 
There  are  major  structural/interface  changes  at  Prof  1.03  and  no  drops  are  detected  with  the  highest  magnification  used 
to  view  these  high  resolution  images.  Thread-  or  finger-like  entities  seen  at  the  interface  are  not  broken  up  into  droplets 
and  are  seemingly  dissolved  at  different  distances  from  the  jet  dark-core  region.  This  forms  a  mixing  layer  in  which 
phase  transition  and/or  large  local  density  noifjuniformities  occur.  A  further  increase  of  chamber  pressure  decreases  the 
length  and  the  thickness  of  the  internal  dark  erne,  and  images  progressively  resemble  injection  of  a  gaseous  turbulent  jet 
into  a  gaseous  environment  (see  magnified  images  in  Fig.  4).  A  gradual  transition  from  a  classical  liquid-like 
appearance  in  which  ligaments  and  drops  are  formed  at  the  interface  in  the  liquid  atomization  regime  to  a  comb-like 
structure  near  the  critical  point,  and  finally  to  where  a  submerged  turbulent  gas  jet  appearance  emerges  can  be 
observed.  Inspecting  a  large  set  of  images  at  high  magnifications,  no  evidence  of  drop  formation  is  seen  in  this  gas-like 
jet  regime.  The  reason  for  this  behavior,  particularly  for  the  change  into  gas-jet  like  behavior,  should  be  sought  in 
progressive  reduction  of  surface  tension  and  heat  of  vaporization  until  they  both  vanish  at  and  above  the  critical  point. 
Similar  observations  are  made  for  the  02-into-N2  case  but  with  the  transition  to  gas-like  jet  at  near  the  oxygen  critical 
pressure  (Pr  =  0.85).  Observation  of  the  gas-like  jet  behavior  and  the  lack  of  any  drop  formation  raises  a  question  on  the 
relevancy  of  current  injection  models  and  some  drop  vaporization/combustion  results  within  this  regime. 


Figure  1.  Schematic  diagram  of  experimental  setup  for  sub-  to  supercritical  jet  injection 


Traditionally,  the  breakup  and  atomization  regimes  of  the  liquid  jet  within  a  narrow  chamber-to-density  ratio  range  is 
shown  on  a  plot  of  Ohnesorge  number  versus  Reynolds  number  (ignoring  gas  density  and  nozzle  geometrical  factors) 
or  by  a  more  refined  criterion  of  Reitz  and  Bracco  [11]  in  which  chamber-to-injectant  density  ratio,  nozzle  geometrical 
factor,  and  Weber  and  Reynolds  numbers  are  all  included.  In  Fig.  3,  as  chamber  pressure  approaches  the  critical 
condition,  surface  tension  is  reduced  to  a  near-zero  value  and  for  the  range  of  Reynolds  number  calculated  here  (25,000 
to  75,000),  the  Ohnesorge  number  sharply  swings  from  a  low  (estimated  to  be  2.8xl0‘3 )  to  a  very  large  (infinity  when 
surface  tension  is  zero)  values.  The  large  Ohnesorge  number  regime  and  its  importance  are  indicated  in  Faeth  [12]  and 
Tseng  et  al.  [13].  At  the  low  Ohnesorge  number  in  this  work,  the  second-wind  induced  liquid  jet  breakup  behavior  is 
observed,  see  Fig.  3.  Past  the  second  wind-induced  behavior,  it  seems  that  before  the  jet  has  the  opportunity  for  full 
atomization,  the  surface  tension  is  sufficiently  and  rapidly  reduced  so  one  achieves  gas-like  jet  appearance  near  but 
before  reaching  the  critical  pressure  value,  see  Fig.  3.  The  aforementioned  threads  or  fingers  are  gasified  and  no  drops 
are  detected.  Transition  into  a  full  atomization  regime  is  therefore  inhibited.  The  investigated  jet  exhibits  liquid-jet  like 
and  gas-jet  like  features  depending  on  the  magnitudes  of  the  surface  tension  and  heat  of  vaporization.  Approximate 
characteristic  time  estimates  by  Chehroudi  et  al  [9]  also  point  towards  inhibition  of  the  atomization. 

JET  GROWTH  RATE 

The  initial  jet  spreading  angle  or  its  growth  rate  is  measured  for  all  acquired  images  and  results  along  with  those  of 
others  are  presented  in  Fig.  5.  The  angles  in  our  work  are  measured  from  the  information  within  a  Mm  distance  close 
to  the  injector  exit  face  (distance-to-diameter  ratio  of  up  to  28).  Chehroudi  et  al  [9]  show  that  the  intact  core  of  the 
liquid  sprays,  Chehroudi  et  al  [  15],  and  the  potential  core  for  different  jets,  Abramovich  [14],  have  comparable  lengths 
to  the  streamwise  extent  of  the  initial  region  used  for  the  growth  rate  measurements,  therefore  ensuring  the  existence  of 
the  mixing  layer.  The  importance  and  justification  of  the  data  sets  selected  in  Fig.  5  and  the  nature  of  their 
measurements  by  other  researchers  are  elaborated  with  sufficient  details  in  our  earlier  paper  to  provide  a  deeper 
appreciation  of  fiw^Fig.  5  and  its  uniqueness,  see  Chehroudi  et  al  [9].  Since  the  jets  investigated  here  exhibit  both 
liquid-like  and  gas-like  jet  appearances,  appropriate  results  for  both  are  presented  in  Fig.  5.  The  simplest  is  the 
prediction  of  the  linear  jet  growth  for  the  turbulent  incompressible  submerged  jet  using  the  mixing  length  concept.  A 
semi-empirical  equation  by  Abramovich  [14]  incorporating  the  effects  of  density  variations  using  a  characteristic 
velocity  is  also  shown,  see  also  Chehroudi  et  al  [9]. 

Brown  and  Roshko  [16]  measure  spreading  for  a  subsonic  two-dimensional  incompressible  turbulent  mixing  layer  in 
which  helium  and  nitrogen  are  used.  They  make  a  distinction  between  mixing  layers  in  which  density  changes  are 
caused  by  temperature  changes,  high-speed  compressible  (high  Mach  number  and  supersonic)  flows,  and  differences  in 
molecular  weights  (i.e^  different  gases).  In  the  jet  we  are  investigating  here  there  are  both  differences  in  temperature  as 
well  as  molecular  weight  except  for  the  N2-into-N2  case.  Their  measurements  are  shown  in  Fig.  5.  Papamoschou  and 
Roshko  [1]  proposed  a  theoretical  equation  for  the  incompressible  variable-density  mixing  layer  as  shown  in  Fig.  5. 
Finally,  Dimotakis  [2]  uses  tire  observation  that,  in  general,  the  entrainment  into  the  mixing  layer  from  each  stream  is 
not  the  same  and,  in  a  system  moving  with  the  aforementioned  convection  velocity,  offers  a  geometrical  argument  to 
derive  an  equation  for  tire  two-dimensional  incompressible  variable-density  mixing  layer,  see  Fig.  5. 

Because  both  liquid-like  and  gas-like  visual  jet  behaviors  are  observed,  the  growth  rate  for  liquid  sprays  produced  from 
single-hole  nozzles  typical  of  the  ones  used  in  diesel  engines  are  also  shown.  Because  of  profound  nozzle  geometrical 
effects,  Reitz  and  Bracco  [11]  and  Hiroyasu  and  Arai  [17],  isothermal-spray  angle  equations  proposed  by  Reitz  and 
Bracco  [11]  at  two  different  nozzle  length-to-diameter  ratios,  along  with  their  corresponding  vertical  error-bands 
indicating  experimental  scatter  around  them,  are  shown  in  Fig.  5.  As  a  cross  check,  a  recent  curve-fitted  equation  to 
experimental  data  from  transient  sprays  proposed  by  Naber  and  Siebers  [18]  is  also  shown.  Their  angle  measurement 
zone  extends  beyond  our  initial  region  and  this  to  some  extent  contributes  to  disagreement  seen  between  the  two  sets  of 
data  for  liquid  sprays  at  injector  length-to-diameter  density  ratio  of  about  4,  see  Fig.  5. 

Figure  5  covers  a  density  ratio  of  three  orders  of  magnitude,  from  liquid  sprays  to  supersonic  mixing  layers,  a  unique 
and  new  plot  on  its  own  right.  Clearly,  within  the  range  plotted,  the, results  of  constant  spreading  angle  of 
incompressible  turbulent  jet  overpredicts  nearly  all  others  in  Fig.  5.  There  also  increasing  disagreement  between 
turbulent  gas  jet  of  Abramovich  [14]  and  incompressible  variable-density  model  of  Papamoschou  and  Roshko  [1]  as 
density  ratio  increases.  To  some  extent,  for  measured  values,  disagreements  in  this  figure  can  be  attributed  to 
differences  in  the  definition  of  the  mixing  layer  thicknesses  and  their  measurement  methods,  see  Chehroudi  et  al  [9].  It 
is  clear  that  for  a  range  of  density  ratios  in  which  our  images  show  gas-jet  like  appearance,  the  experimental  data  agrees 
well  with  the  proposed  theoretical  equation  by  Papamoschou  and  Roshko  [1]  and  closely  follows  the  trend  of  Dimotakis 
[2].  This  can  be  taken  as  further  and  quantitative  confirmation  that  at  ambient  supercritical  pressure  and  temperature 
conditions  (based  on  the  pure  injectant),  the  injected  jets  visually  behave  like  a  gas  though  technically  it  may  be 
referred  to  as  “fluid”.  To  our  knowledge,  this  is  the  first  time  such  quantitative  verification  has  been  demonstrated. 


Above  the  critical  condition,  there  is  marked  disagreement  in  both  magnitude  and  slope  between  liquid  sprays  (at  a 
comparable  length-to-diameter  ratio  of  85)  and  our  data,  see  Fig.  5,  even  though  the  jet  investigated  here  appears  to  go 
through  initial  phases  of  the  liquid  atomization  process,  see  Figs.  3.  The  reason  is  that  although  the  jet  studied  here 
shows  second  wind-induced  breakup  features  similar  to  liquid  jets,  it  fails  to  reach  full  atomization  state  as  chamber 
pressure  (really  density)  is  raised.  This  is  because  the  thermodynamic  state  approaches  the  critical  point  and 
consequently  both  surface  tension  and  heat  of  vaporization  are  reduced  to  near-zero  values.  Transition  into  the  full 
liquid  atomization  region  is  therefore  inhibited. 

FRACTAL  ANALYSIS 

In  the  past  4b  years,  a  number  of  applications  of  fractal  analysis  have  been  demonstrated  in  different  disciplines.  For 
example,  Sreenivasan  and  Meneveau  [19]  compute  the  fractal  dimension  of  the  turbulent/non-tufbulent  boundary  of  an 
incompressible  axisymmetric  gaseous  jet  and  report  a  value  of  1.33  using  two-dimensional  slicing  by  a  laser  sheet. 
Sreenivasan  [20]  also  reports  values  of  1.35,  1.34,  and  1.38  for  a  round  gaseous  jet,  a  plane  gaseous  mixing  layer,  and  a 
boundary  layer  flow,  respectively.  Figure  6  shows  results  of  plotting  the  fractal  dimension  applied  to  the  N2-into-N2  jets 
as  a  function  of  relative  chamber  pressure.  In  this  figure  the  fractal  dimension  by  other  researchers  in  the  gaseous  jets, 
mixing  layer,  and  boundary  layers  are  also  shown  for  comparison  purposes.  The  interesting  part  is  diat  the  average 
value  measured  for  our  jet  reaches  very  near  those  of  gaseous  jets  and  mixing  layers  above  the  critical  point  of  the 
injectant.  This  is  additional  quantitative  confirmation  of  the  gaseous  jet  behavior  indicated  through  our  growth  rate 
measurements  in  this  region.  As  pressure  (density)  is  decreased  below  Pr  of  about  0.6,  the  fractal  dimension  is  rapidly 
reduced  towards  the  Euclidean  value  of  1  for  a  smooth  circular  cylinder  with  no  surface  irregularities.  Considering  the 
geometric  interpretation  of  the  fractal  dimension,  measured  values  are  quite  consistent  with  the  visual  observations  of 
jet  interface  changes  in  our  tests.  For  more  details  see  Chehroudi  et  al.  [21]. 

SUMMARY  AND  CONCLUSIONS 

Anatomical  changes  and  growth  rates  of  jets  injected  into  an  environment  at  fixed  supercritical  temperature  but  varying 
pressure  from  sub-  to  supercritical  are  analyzed.  As  chamber  pressure  is  increased  from  a  low  subcritical  value,  the  fluid 
in  the  jet  appears  to  go  through  classical  liquid  jet  breakup  stages  up  to  a  second  wind-induced  breakup  regime  with 
ligaments  and  many  droplets  ejecting  from  tire  jet.  Penetration  into  the  full  atomization  regime  is  inhibited  near  but 
before  the  critical  pressure  of  the  injectant  because  of  the  combined  effects  of  lowered  surface  tension  and  heat  of 
vaporization.  At  this  point  the  jet  assumes  a  gas-jet  like  appearance  that  remains  up  to  the  highest  pressure  tested  here. 
Also,  a  unique  and  comprehensive  plot  on  growth  rate  is  formed  using  most  relevant  works  of  others  covering  an 
ambient-to-injectant  density  ratio  range  of  three  orders  of  magnitude.  Our  measured  jet  growth  rates  follow  theoretical 
equations  proposed  by  Papamoschou  and  Roshko  [1]  and  Dimitakis  [2],  both  derived  for  incompressible  variable- 
density  gaseous  turbulent  mixing  layers,  starting  at  a  pressure  near  but  below  the  thermodynamic  critical  pressure  of  the 
injectant  quantitatively  supporting  the  observed  gas-jet  like  visual  appearance  of  the  supercritical  jets  for  the  first  time. 
This  jet  exhibits  both  liquid-jet  like  and  gas-jet  like  faces  depending  on  the  values  of  surface  tension  and  heat  of 
vaporization.  The  appropriateness  of  the  surface  tension  for  liquid  jets  and  sprays  and  its  irrelevance  for  gaseous  jets  are 
among  the  issues  to  be  reconciled.  The  first  fractal  analysis  of  the  jet  for  the  supercritical  case  also  shows  agreement 
with  values  reported  for  gaseous  jets/mixing  layers,  another  quantitative  indication  of  the  gas-jet  like  behavior. 
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Figure  3.  Back-illuminated  images  of  nitrogen  injected  into  nitrogen  at  a  fixed  supercritical  temperature  of  300  K  but 
varying  sub-  to  supercritical  pressures  (Pcnucai  -  3.39  MPa;  Tc  =  126.2  K).  Pch/Pmucai=  2.74,  2.44,  2.03,  1.62,  1.22,  1.03, 
0.83,  0.62,  0.43,  0.23;  from  upper  left  to  lower  right.  Re=  25,000  to  75,000.  Injection  velocity:  10-15  m/s.  Fraud 
number:  40,000  to  1 10,000  (momentum-dominated  jet).  Injectant  temperature:  99  to  1 10  K. 
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Figure  4.  Software  magnified  images  of  the  jet  at  its  outer  boundary'  showing  transition  to  the  gas-jet  like  appearance 
starting  at  just  below  the  critical  pressure  of  the  injectant.  Images  are  at  fixed  supercritical  chamber  temperature  of  300 
K. 
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Figure  5,  Shows  spreading  or  growth  rate  as  tangent  of  the  visual  spreading  angle  versus  the  chamber-to-injectant  density  ratio.  (*) 
refers  to  data  taken  at  AFRL. 
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Figure  6.  Box-counting  and  Minkowski  (Euclidean  Distance  Mapping  (EDM)  algorithm)  fractal  dimensions  of  the 
visual  boundary  of  the  jet  as  a  function  of  the  reduced  chamber  pressure  for  N2-into-N2  injection. 


